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ABSTRACT 

We present a new framework for the extraction of the strong couphng from 
hadronic r decays through finite-energy sum rules. Our focus is on the small, but 
still significant non-perturbative effects that, in principle, affect both the central 
value and the systematic error. We employ a quantitative model in order to 
accommodate violations of quark-hadron duality, and enforce a consistent treat- 
ment of the higher-dimensional contributions of the Operator Product Expansion 
to our sum rules. Using 1998 OPAL data for the non-strange isovector vector and 
axial-vector spectral functions, we find the rif = 3 values as(m^) = 0.307±0.019 
in fixed-order perturbation theory, and 0.322±0.026 in contour-improved pertur- 
bation theory. For comparison, the original OPAL analysis of the same data led 
to the values 0.324 ± 0.014 (fixed-order) and 0.348 ± 0.021 (contour-improved). 



* Permanent address: Departament de Fisica, Universitat Autonoma de Barcelona, E-08193 Bellaterra, 
Barcelona, Spain 
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I. INTRODUCTION 



In the past few years there has been a renewed interest in the precision determination of 
tts from non-strange hadronic r decays. One reason for this interest is the recent calculation 
of the coefficient of the O(a^) term in the perturbative contribution to the Adler function 
[l|. This contribution dominates the ratio of the hadronic r decay width and the electronic 
decay rate pi], 

Rr = r[r^ — z/,- hadrons (7)]/r[r^ Ure'Veil)] ■ (1-1) 

Another reason is the existence of a number of competing analysis methods which lead to 
results that are not, or only barely, consistent with one another. In fact, the error on as 
from r decays quoted in a recent (2009) review p.] has gone up since its 2006 version, for the 
simple reason that the result of Ref. p] was obtained by averaging the central values of all 
recent r decay determinations of and the error by considering the spread of these central 
values. All determinations are based on data from (primarily) the ALEPH (see Ref. j3] for 
their 1998 analysis and Ref. jsj for their 2005 analysis) and (also) the OPAL (see Ref. 0) 
collaborations; the differences in the results are on the theory side. 

Clearly, this is an unsatisfactory situation. There are at least three theoretical issues 
related to the discrepancies between the different determinations. Number one is the long- 
standing question as to which resummation scheme, fixed-order perturbation theory (FOPT) 
or contour- improved perturbation theory (CIPT), is best used for evaluating the perturbative 
contributions to Rr- Many of the recent reanalyses have focussed on this question 0- 



Much less attention has been devoted to the two other issues, both of which concern 
non-perturbative contributions to R^-- Because of the relatively low value of the r mass, 
such contributions cannot be entirely neglected, even if they are expected to be small. Issue 
number two is the question of whether the Operator Product Expansion (OPE) contributions 
beyond perturbation theory have been consistently taken into account. Here we are aware 



of only one systematic investigation [12|, in which it was demonstrated that self-consistency 
problems existed for a number of earlier analyses. Specifically, it was shown that the OPE 
parameters obtained in those analyses do not provide a good match between data and theory 
when the upper limit Sq on the hadronic invariant mass-squared, s, in the weighted integrals 
over the spectral functions (which enter the finite-energy sum rules (FESRs) employed in 
these analyses) is varied away from m^} Issue number three concerns potential violations of 
"quark-hadron duality" not taken into account in previous FESR determinations of a^. In 
the case of FESRs, the assumption of quark-hadron duality amounts to presuming that all 
non-perturbative effects are accounted for by higher-dimensional terms in the OPE. To date. 



this assumption has not seen a systematic investigation. As already explained in Ref. [12| , the 
issues of the OPE and possible violations of quark-hadron duality are intricately connected: 
without a quantitative analysis of duality violations, it turns out to be difficult to treat 
the OPE consistently without relying on "external" results. For instance, in the analysis 



of Ref. [12[ , the dimension-4 term in the OPE had to be fixed using a result for the gluon 
condensate from charmonium sum rules. 

In this article, we aim to address this situation, focussing on the non-perturbative ques- 



Some of the earlier analyses carried out this test for the FESR based on the kinematic weight (the weight 
yielding when sq — m"^), for which it works reasonably well. Reference jl^l showed that these analyses 
unambiguously fail the tests for other doubly-pinched FESRs. 
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tions. We use a recently developed model for the duality- violating (DV) part of the wd-flavor 
vector (V^) and axial- vector (A) spectral functions [l3|, which makes it possible to carry out 
a self-contained FESR analysis in which stability with respect to varying sq is checked self- 
consistently without reliance on external values for any of the OPE parameters. 

Since no QCD-based theory of quark-hadron duality exists, we have to resort to a model. 
This means that our results will be based on the (testable) assumption that this model gives 
a good description of the DV part of the spectral functions for values of s from s — )■ oo down 
to a minimum value, Smin, sufficiently low to lie in the range s < kinematically accessible 
in hadronic r decays. We emphasize that the need to make such an assumption has always 
been a fundamental "shortcoming" of the determination of from r decays. Assuming 
quark-hadron duality a priori, and therefore neglecting the effect of DVs altogether, also 
amounts to employing a - probably worse - model. In other words, in order to investigate 
the systematics related to the assumption of quark-hadron duality, one cannot avoid the 
adoption of a model of the DV part of the spectral functions. In Ref. [l3| it was found that 
the model we intend to employ gives a reasonable description of the spectral functions in 
the region 1.1 GeV^ < s < rn^? Moreover, the physics of our model is based on a picture of 
the hadron resonances which are experimentally seen in the spectral function. Resonances 
are not described by perturbation theory or the OPE, and thus should be part of any model 
aiming to describe violations of quark-hadron duality. 

In this work, we do not address the issue of the optimal choice of resummation for the 
truncated perturbative series in a given FESR. While this systematic, of course, forms a 
potentially important part of the final theory error on a^, we have no new elements to 
add to the discussion of this issue. Moreover, we believe that the non-perturbative part of 
the systematics should be understood first, in order to get a more reliable picture of the 
quantitative discrepancy between results based on FOPT and CIPT. We will therefore carry 
out our whole analysis with both resummation schemes. 

To date, two experiments, ALEPH 0, Isj and OPAL [6|, have made the non-strange V 
and A spectral functions from their r-decay analyses publicly available. The 2005 analysis 
of ALEPH is more recent, and based on more statistics, and thus would be expected to 
have smaller experimental errors. Unfortunately, the 2005 ALEPH data cannot be used at 
present, because correlations due to unfolding have been inadvertently omitted in the original 



ALEPH analysis and hence from the publicly posted covariance matrices [IJ]. Since the 
re-analysis of the ALEPH data has yet to be completed, we limit ourselves, in this article, 
to an analysis employing the OPAL data. 

In order to normalize the various exclusive-mode components of the spectral functions, 
OPAL relied on the branching fractions available in 1998, as well as the then-current values of 
Vud and the electronic branching fraction B^. All of these have been updated since then, and 
this makes it possible to at least partially update the OPAL inclusive spectral distributions 
as well. Such an update would allow an updated, though still OPAL-based, determination 
of Us- Here, since our primary goal is to investigate the impact of the novel features of 
our treatment of non-perturbative effects on the extracted results for a^, we choose not to 
perform this update, and instead work with the data in precisely the same form as used by 



OPAL [6[. We plan to devote a separate article to an adaptation of OPAL data to recent 



^ In the present, more detailed analysis, we will find that a significantly larger value of s,„i„ is preferred. 
^ We thank members of the ALEPH collaboration for private communications, in which the existence of 
this problem has been confirmed. 
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values of the exclusive branching fractions, Vud and Be, and an investigation of the effect of 
this adaptation on the value of ag and other OPE parameters. 

This article is organized as follows. In Sec. [TT]we present a brief review of the application 
of FESRs to hadronic r decays, with emphasis on the issue of quark-hadron duality. In 
Sec. mil we are then able to provide a more thorough discussion of the various systematic 
errors discussed already above. Preparing for a presentation of our results in Sec. IVTI we 
describe the theory parametrization we will employ in more detail in Sec. \IV\ and discuss 
the issue of strong correlations in the integrated data, and our resulting fitting strategies, 
in Sec. |Vl Apart from reporting on our fits in Sees. IVI Al and IVI Bt we consider also, in 
Sec. IVI Dj the V + A channel sum (related to the non-strange part of Rr) as well as the 
V — A channel difference. In the latter case, we demonstrate that our fit results satisfy the 
Weinberg sum rules 15| as well as the DGMLY sum rule for the vr electromagnetic mass 



difference 16|. Section [VIII contains a summary of our results, including a conversion of 



to its value at the Z mass; Sec. IVIIII contains our conclusions. 



II. THEORY SUMMARY 

Our analysis will involve the correlation functions 

= ^ I rf'xe*^^(0|T{j^(x)jt(0)} |0) (2.1) 

where is one of the non-strange V or A currents, U'f^d or U7^75(i, and the superscripts 
(0) and (1) label the spin. The decomposition in the third line employs the combinations 
n(i+o)(g2^ and g^n^°^(g^), which are free of kinematic singularities. Defining s = = —Q^ 
and the spectral functions 

p(i+o)(5) = il„,n(i+o)(^) ^ (2.2) 

TT 

Cauchy's theorem and the analytical properties of 11*^^+°^ (s), applied to the contour in Fig. [H 
imply the FESR relation 

I^-\{s,) ^ - r dsw{s) p^^;^\s) = --^ / dsw{s) U^^;^\s) , (2.3) 

Sq JO ZlilSQ J\s\=so 

valid for any sq > and any weight w{s) analytic in the region of the contour p/7|]. In the 
present work we will restrict ourselves to polynomial weights. Partial integration allows the 
right-hand side of Eq. (12. 3 p to be recast in terms of the Adler function 

D{s) = -s ^f^i_i£l . (2.4) 

ds 

The spectral functions PyJa\s) are measurable in hadronic r decays. Explicitly, for 
Standard Model decays induced by the flavor ud (isovector) currents, with wt{s; Sq) = 
(1 — s/so)^(l + 2s/so), wl{s; Sq) = 2(s/so)(l — s/soY, and the scaled, non-strange V and A 
widths 

D _T[t- ^ z/^ hadronsy/A;w (7)] 

^V/A;ud = T^f^ ^ / XI ; l^-^j 
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FIG. 1: Analytic structure of Il^^^^\q'^) in the complex s = plane. The solid curve shows the 
contour used in Eg. i2.3]) . 



one has 18 



R 



V/A;ud 



Rv/A;udiso = ml) 



(2.6) 



in which R 



V/A;ud 



[sq) is defined by 



R 



V/A;ud 



So) = l27r'^\Vud\'^SEw — 



So Jo 



so 



ds 



WT{s]So)pyJ^\s) -WL{s]So)py'/A{s)\ , (2.7) 



(0) 



where Sew is a short- distance electroweak correction. Since, in the Standard Model, 



p'^\s) 



0[{md - m„)^)] , 
pf{s) = 2f:{6{s-ml,)-6is))+0[im, + m^)'] 



(2.8) 



the differential distributions proportional to the expression in square brackets in Eq. (12. 7p 
provide a direct measure of pyj'^\s), up to numerically negligible O (m^^^ corrections.'' 
The second delta function in Eq. (12. 8p . which comes from the kinematic singularity present 
in does not contribute to the integral in Eq. (12. 7p as a result of the factor of s in the 
accompanying weight wl(s; Sq). 

For sufficiently large sq, and ignoring for a moment that the OPE is only valid for large 
euclidean Q^, the right-hand side of Eq. (12.30 can be approximated using the OPE for 
s). Experimental spectral data can then be used to fit the OPE, and extract pa- 



n 



(1+0), 
V/A 



rameters such as In what follows, we will denote the experimental version of the 



The TT decay constant, = 92.21(14) MeV 19], is presently known very accurately. The central value of 
the T — > TTj^T- branching fraction, _B^, employed in the 1998 OPAL analysis corresponds to the somewhat 
larger value 94.0 MeV. In order to match exactly the OPAL treatment of the V and A spectral functions 
we employ the latter value in our analysis. Note that, since B^^ was obtained by the PDG in a combined fit 
to the full set of r basis modes, it would, in fact, be inconsistent to change just -Btt without simultaneously 
changing all other branching fractions. We will revert to the updated value /^r = 92.21(14) MeV in our 
later analysis employing updated OPAL data. 
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spectral integral on the left-hand side of Eq. (12. 3 p by Iv/A-exi^o) (generically, /^Jf^(so)) and 
the theoretical representation of the contour integral on the right-hand side by ly/A-thi^o) 
(generically, Ith\so))- 



In the upper part of the energy region allowed by r-decay kinematics 



tt(1+0). 



IS 

OPE 

dominated by its dimension D = contribution, i.e., the perturbative contribution in the 
chiral limit. ^ The perturbative expression for the Adler-function (12. 4p . which is known to 
order can be written as^ 

1 oo n+1 / \ ^=-1 

P(^)]g™ = I-^E <if^')T. kcn J log — ] , (2.9) 

n=0 k=l \ / 

with as(/i^) = as(/i^)/7r. Since D{s) is independent of /i^, we can choose /i^ = —s, indicating 
that only the coefficients c^i are independent; all other Cnk can be expressed in terms of the 
Cni through the renormalization group. In the MS scheme, cqi = cn = 1, Coi = 1.63982, 
C31 = 6.37101 and C41 = 49.07570 [1]. We will use the guess C51 = 283 of Ref. Q for the next 
coefficient, assigning an uncertainty of ±283 in order to estimate the error due to truncating 
perturbation theory (c/. Sec. IVI C\i . 

The freedom to choose /x^ in Eq. (12. 9 p is at the heart of the different prescriptions employed 
for evaluating the perturbative contribution to the right-hand side of Eq. (12. 3p : in FOPT, 
/i^ = So is used in Eq. (12. 9p . whereas in CIPT /x^ = —s is employed inside the contour 



integral on the right-hand side of Eq. (12. 3 p |21 



Beyond perturbation theory, one may improve the approximation to the right-hand side 
of Eq. (12. 3p by including higher-dimension contributions to nopE'*(s). Explicitly, 

nSE°'W = f;^, (2.10) 

k=0 V ^) 

with the OPE coefficients C2k logarithmically dependent on s through perturbative cor- 
rections. The term with k = corresponds to the purely perturbative, mass-independent 
contributions, represented by Eq. (12. 9p . The C2k, with k > 1, contain non-perturbative 
D = 2k condensate contributions, and are, in principle, different for the V and A channels. 

We will neglect C2, which is purely perturbative and quadratic in the light quark masses.^ 
It has been suggested that a non-perturbative D = 2 term should be added in order to 
account for the truncation of the perturbative series for the D = term [22|. We postpone 
an investigation of this issue to future work, and here set C2 = 0. 

Neglecting contributions of 0{m'l^^ or proportional to (uu) — (dd), both of which are 
numerically very small, the coefficient C4 is a linear combination of the "gluon condensate" 
{asG^yG^") (with G^u the gluon field strength), and the chiral condensates mi{ip^ipi), i = 
u, d, s. To leading order in as both contributions to C4 are the same in the V and A 
channels. Differences in the D = 4 V and A light quark condensate contributions enter 
beginning at 0{as) or 0{m'l^^ [2^ . 



^ Perturbative contributiaiis proportional to powers of the quark masses are included in Z? > 2 OPE terms. 
^ See, for instance, Ref. [20|. 

^ Since in the present study we are only dealing with the light up- and down-quark correlators, the D = 2 
mass-squared corrections are tiny. Still, one version of our analysis code has implemented all known 
corrections up to 0{al), and we have verified that they are indeed negligible. 
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The coefficient Cg is assumed to be dominated by four-quark condensates, because the 
contribution from [g^ fabcG'^^G^^G'^^) vanishes at leading order in as (2^1, and the contribu- 
tions from lower-dimensional operators are suppressed by powers of the quark mass. Coef- 
ficient functions for the four-quark condensates were calculated to next-to-leading order in 
Ref. M. Ki D = d> there is a proliferation of operators, and very little detailed information 



is available. As we will explain in Sec. IIVI below, we will not need to consider terms with 
L> > 8. 

The OPE is valid when the euclidean distance |x| in Eq. (12. ip is small compared to ^^CD^ 
or, equivalently, when euclidean is positive and large. However, both perturbation theory 
and the OPE are expected to break down near the positive real axis |26| . We may account 
for this additional non-perturbative effect by writing the right-hand side of Eq. (12.31) as [l3[ 

with 

A(.) = n(i+°)(.)-ng;°)(.). (2.12) 

The difference A(s) defines the duality violating contribution to 11'^^"''°) (s). 

All previous determinations of from hadronic r decays have assumed, implicitly or 
explicitly, that integrated DV contributions are small enough to be neglected for the weights 
employed in the analysis. While this assumption has sometimes been checked for self- 



consistency (see, e.g., Ref. [12|), a comprehensive quantitative estimate of the impact of 
DVs on the precision with which can be determined has not been provided. One of 
the aims of the present work is to provide a comprehensive analysis which takes DVs into 
account, and hence remedies this shortcoming. 



As shown in Ref. [13|, if A(s) is assumed to decay fast enough as s — )• oo, the right-hand 



with 



side of the FESR relation (12. 3p can be rewritten as 

ll:\so) = / dsw{.s)n^^tl\s)+V^{so) , (2.13) 

Vu,iso) = / dsw{s)-lmA{s) . (2.14) 

Sq JsQ TT 

The imaginary parts ^ lmAv/A{s) can be interpreted as the DV parts, Py/^(s), of the V/A 
spectral functions. Following Ref. |l3[, we will parametrize PvJaI^) 

P°7a(s) = ^v/A e-^^/-^^ sin {ay/ a + I^v/as) ■ (2.15) 

This introduces, in addition to and the D > 4 OPE condensates, four new parameters 
in each channel. The OPE, supplemented by the ansatz f l2.15p . will be assumed to hold 
for s > Smini where Smin will have to be inferred from fits to the experimental data. The 
extended analysis, including DVs, will of course only be possible if Smin lies significantly 
below m^. 

The ansatz (I2.15P is modeled on the asymptotic behavior for large s of a semi-realistic 



model for the QCD spectrum in a given channel. This model was developed in Ref. |27 
based on earlier ideas described in Ref. j28|. It incorporates a combination of large- iVc 
insights (narrow resonances with widths increasing with mass) as well as the Regge picture 
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for the spacing between resonances. These ingredients lead quite naturally to the exponential 
decay in Eq. f l2.15p with the decay parameter 7 ~ 1/iVc, as well as the oscillatory behavior 
represented by the sine function, both with arguments (approximately) linear in s. We 
favor this model over other attempts to model DVs because of its natural connection to 
the resonance structure of the spectral distributions, something that is not evident in other 
model s ( such as those based on instantons). For detailed discussions of the model 



see 



Refs. |27|, l29|, |30 



III. SYSTEMATIC ERRORS 

There are three sources of systematic error affecting, to various extents, existing FESR 
determinations of a^. Since the investigation of two of these sources is the central aim of 
this work, we briefly describe each of the three sources here, before embarking on the details 
of our analysis. 

1. There are (at least) two ways of partially resumming the perturbative contribution 
to lj:^\sf)), CIPT and FOPT (c/. Sec. [TTl)- The relative merits of the two methods 
have been the subject of a number of investigations [l|, 0"11|. While no particular 
preference is given to either scheme in Refs. [H, llH, CIPT is favored in 0, [sl, whereas 



Refs. [8|, Il0| give arguments in favor of FOPT, in the latter work through a new CI 
expansion in a conformally mapped coupling. 

We will not attempt to resolve the associated systematic uncertainty in this work, 
but instead report on the results of our fits using both CIPT and FOPT. In fact, it 
is interesting to see what discrepancy remains between CIPT and FOPT after other 
systematic errors, described below, have been properly taken into account. 



2. With the exception of Ref. [12j, the OPE has not been treated consistently in previous 
extractions of as from r decays, in the sense we now explain. Consider a term of 
order in the OPE of Eq. fl2.1UI) . The dominant term in the expansion in 
of the corresponding coefficient Cd, D = 2k, is a constant of order one (times the 
relevant condensate). In the sum rule (12.31) . this term in the OPE is picked out by 
the term of degree /c — 1 in the weight w{s). Other terms in w{s) will also pick up 
contributions from Cd{s), because of the logarithmic dependence of Cd{s) on s, but 
such contributions will not be dominant as they are suppressed by at least one extra 
power of as- Thus, if weight functions up to degree n are used in the fits, it follows 
that terms up to at least order k = n + 1 must be kept in the OPE in order to retain 
all potentially relevant contributions not suppressed by at least one extra power of a^. 

In the conventional analysis of Refs. 0-0, Sff this was not done: while weights up to 



degree n = 7 in s were employed (which would generally require keeping terms in the 
OPE up to D = 16 (/c = 8)), only terms up to D = 8 (/c = 4) were retained. As noted 
earlier, it was found in Ref. that this is not self-consistent: with the parameter 
values found in those fits, the sq dependence of the theory curves does not match that 
of the data for the majority of the weights employed, as well as for alternate degree 2 
and 3 weights which explicitly test the D < 8 parameters obtained from the original 
fits.^ In this work, we will restrict ourselves to weight functions of degree n < 3, 



See also Ref. 



27[, where it was shown in a model study that adding a. D — 10 term to the OPE in an 
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corresponding to keeping terms up to D = 8 in the OPE. As we will explain below, 
this implies we will have to vary sq in the FESR (12.31) : it is not possible to restrict a 
consistent analysis to only sq = m"^ without additional uncontrollable assumptions. 

3. Typically, previous analyses neglected the presence of duality violations.^ While in 
some cases this assumption was checked for self-consistency |12|, such a check does 
not provide a quantitative assessment of the impact of residual DV effects. It has long 
been known that FESRs with a simple weight like w = 1 have sizable DVs, even at 
scales ~ 2 — 3 GeV^, but that switching to weights "pinched" (having a zero) near 
s = Sq significantly reduces this effect 3^. It has become a standard assumption that 



using only weights which are at least doubly pinched will suppress DVs sufficiently so 
as not to affect the value of, or error on, a^. Clearly, in view of the quite small errors 
on as reported in the recent literature, this issue is in need of further investigation. 

The use of weights which are at least doubly pinched, i.e., which contain at least two 
powers of Sq — s, forces us to vary sq in Eq. (12. 3p . if the OPE is to be treated consistently 
in the sense described above. Suppose one wants to consider only sq = m^, and fit as 
as well as C4^,„^D=2k, i-e., k parameters. To fit k parameters, one requires at least 
k + 1 data points, and therefore at least k + 1 linearly independent weights if one 
uses only sq = m^. If all weights contain the factor (sq — s)^, the minimally required 
highest degree will he n = k + 2. This is inconsistent with our criterion of point 2 
above, which would require terms up to order + 3 to be kept in the OPE. The only 
way out is to vary sq, and/or to consider weights that are less than doubly-pinched. 
This makes the need to take DVs into account more urgent, because it is certainly not 
justified to assume that integrated DVs are negligible for weights which are less than 
doubly pinched, over any sizable interval in sq below m^. 



IV. PARAMETRIZATION USED IN FITS 



In this section, we describe in detail the parametrization of the theory that we will use 
in our analysis of the data. 



A. Selection of moments 

We wish to consider terms in the OPE only up to D = 8, for several reasons. First, we 
expect that small contributions to typical OPE integrals associated with these condensates 
will be potentially sensitive to residual integrated DV contributions, making it possible to 
check the impact of DVs on earlier determinations of the condensates. Second, it appears 
unlikely that we can reliably determine condensates with D > 8 from existing data. Hence 
we focus on FESRs where such contributions will be strongly suppressed. Finally, little is 
known about the OPE, but it is almost certainly not a convergent expansion. One would 



analysis like that of Refa-JpriZl can make a significant difference. 
^ Exceptions are Refs. [7|, |32|; however, as is clear from the results of the present work, those investigations 

of DVs involved additional assumptions which we are able to avoid in the present analysis. 
° Recall that we set C2 = in this article. 
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thus expect it to break down at some sufficiently high order, making it prudent to hmit 
ourselves to a relatively low maximum order. 

From the arguments in Sec. Illlt this restricts us to weights with degree < 3, i.e., to at 
most four linearly independent weights. Since we will be fitting up to four OPE parameters 
{as as well as the D = 4, 6 and 8 condensates), in addition to the DV parameters in each 
channel, this already forces us to consider the sum rules fl2.3l) with more than one value of Sq 
for at least one weight. We will vary sq over the interval [smin, f^l], and explore the stability 
of the fits as a function of Smin- 

In this work, we choose to consider the weights^^ 

M^) = 1 , (4.1) 

W2{x) = 1 — , 

wsix) = {I - xY{l + 2x) = Wt{s] So) , 

X = s/sq . 

A key point is that we explicitly incorporate DVs in our fits, and therefore need to use at least 
one weight sensitive not only to Og and the OPE coefficients, but also to the DV parameters. 
This stands in contrast to other work to date, where a desire to neglect DVs motivated the 
use of (at least) doubly-pinched weights, which are known to suppress such contributions. 
Such doubly-pinched weights are, in fact, too insensitive to the DV parameters to allow for 
reliable ffis of these parameters. In order to maximize our sensitivity to DVs and hence 
improve our ability to fit DV parameters we include the unpinched weight Wq in all our fits. 
The other weights have been chosen by requiring them to be of degree < 3, to have no term 
linear in s, and to be singly pinched {1112) or doubly pinched (t&s). Alternative sets of weights 
satisfying the same requirements are obtained by replacing either W2 or with 1 — x^. Our 
results with these alternative sets are completely consistent with those obtained from the 
set dUD. 

For constant Co, Ith^"^ picks out the D = 6 term, while /j^"^'' picks out the D = 6 and 
D = 8 terms. In practice, the logarithmic dependence of Co on s beyond leading order in 
as implies that all terms in the OPE contribute for all choices of w{s). However, Cq and 
Cs will be primarily determined by the W2 and W3 FESRs. In the present work, we will 
represent the D > 4 contributions using effective values C4, Cq and Cg independent of s. 
This implies that C4 does not contribute to ffis involving any of the moments of Eq. fl4.1l) . 
For the case of C4, we have checked, in the case of ffis to /^"(sq), that the numerical effect 
of this approximation is tiny, cf. Sec. IVI A[ 

As we will show in Sec. IVTl ffis to moments with weights Wq, W2 and W3 lead to stable and 
self-consistent results. One might also consider including moments with weights containing 
a linear term in s, such as the weight wi{x) = 1 — x. Such moments are sensitive to the gluon 
condensate and would thus allow us to estimate its size, although due care would have to be 
exercised with the interpretation of such a result since this condensate mixes with the unit 
operator. However, renormalon- inspired model studies of higher orders in perturbation 
theory along the lines of Refs. [sl, 20] appear to indicate that perturbation theory, truncated 
at currently known orders, may be less well converged to the full resummed result for iI^^\sq) 
than for /^^^"'^'^^(sq). This may be related to the observation that the = 4 term, present in 



We use hats in order to distinguish our set of weights from a different set of weights considered in Ref. 12|. 
A similar observation holds for some of the condensates in Cq and Cg. 
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the OPE representation of Jji^ (sq), is affected by tlie leading infrared-renormalon ambiguity 
in tlie Borel resummation of the associated perturbative series. In contrast, non-perturbative 
contributions to the OPE representations of I^^^'^\so) depend most significantly on the Co 
with D > 4, and are affected primarily only by subleading ambiguities, associated with 
more distant infrared renormalon poles. This is one reason we restrict ourselves to the 
weights fl4.ip in this article. 

We have also studied weights with a term linear in x, such as Wi{x), added to our set of 
weights (14. ip . but find that, with analyses of the type we present in Sec. IVll the uncertainty 
on the resulting determination of the gluon condensate remains large compared to that of 
other determinations in the literature. We intend to further investigate such fits, and in 
particular the determination of the gluon condensate, in a forthcoming article. Here we just 
note that fits similar to those presented in Sec. I VI I but also including the weight Wi yield 
results consistent with the fits presented in Sees IVI Al and IVIBI Finally, we observe that 
■"^0,2,3 (a;) are a complete, linearly-independent set of weights of degree < 3 without a term 
linear in s. 

B. Duality violations 

We will parametrize the duality- violating part of pv and pa as in Eq. f l2.15p . This 
introduces four new parameters in each channel, forcing us to consider values of Sq over an 
interval [sminy^r]- Since our ansatz f l2.15p is assumed to hold only for sufficiently large s, 
Smin must be large enough to lie in the region of validity of this assumption, but low enough 
to be kinematically accessible in hadronic r decay. Our ansatz is therefore only practical if 
we also assume that such an Smin exists. We are then interested in a value of Smm that is 
small enough to maximize the data available for use in our fits, requiring at the same time 
that the DV ansatz with that choice of s^m gives a good description of the data. 

A priori, there is no reason for s^m to be equal in the V and A channels. We therefore 
present two types of analysis: one for the V channel alone, and one for the combined V and 
A channels {VhA). In the latter case, we will employ an Smin common to both channels. 
This is equivalent to assuming an Smin < ""^r exists such that the asymptotic behavior has set 
in for both the V and A channels for all s > Smin- Since it seems unlikely for the asymptotic 
behavior in a given channel to set in below the lowest resonance in that channel, we expect 
to find Smin ~ '^ai higher for the combined VhA fits. In practice, we find an optimal 
choice Smin ~ 1.4 — 1.5 GeV^. 

We have also considered fits to only the A channel, but find that the data in that channel 
lead to a poor determination of the DV parameters, and thus also of Us- We believe this 
is due to the lower quality of the A-channel data, rather than the absence of a sufficiently 
low A-channel Smm, but it is impossible to decide this from the data alone. We therefore do 
not discuss purely A-channel fits, and restrict our analysis of this channel to combined fits 
involving the V channel as well. 

V. CORRELATIONS AND FITTING STRATEGIES 

Values of the left-hand side, J^Jf^(so), of Eq. (12.30 for nearby values of sq are very strongly 
correlated, and this has repercussions for the choice of fitting strategies. We describe the 
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strong correlations in Sec. IV At and our strategies in Sec. IV Bl below. As already explained 
in the introduction, we limit ourselves to an analysis of the OPAL data j6|. 

A. Correlations and errors 

The data we will use are the OPAL compilation of the non-strange V and A spectral 
functions. ^'^ These data appear in our fits through the weighted integrals, /^Jf^(so), appearing 
on the left-hand sides of the FESRs (12.31) . for the various weight functions we consider. 
These integrals are, of course, represented numerically by sums over the appropriate sets of 
experimental bins. Since OPAL's bin width is 0.032 GeV^, varying sq between approximately 
1.5 GeV^ and 3.120 GeV^ (OPAL's highest bin in the V channel, and almost equal to 
ml = 3.157 GeV^) or 3.088 GeV^ (OPAL's highest bin in the A channel) provides about 50 
data points for each integral. 

The integrals /^Jf^(so) are, however, highly correlated. For instance, if we consider /^°(so) 
on the interval sq G [1.504,3.136] GeV^,^^ the corresponding correlation matrix has sub- or 
super-diagonal elements as large as 0.998, a largest eigenvalue ~ 33, and a smallest eigenvalue 
~ 0.00019. Nonetheless, as we will show in Sec. IVI At it turns out that reliable, standard 

fits to the data for wq, using our parametrization of the right-hand side of Eq. (12.31) . are 
possible. 

The situation changes if we consider two or more weight functions simultaneously. Fo- 
cussing on our primary fits, with moments constructed with the weight functions i&o,2,3; not 
only the correlations for each moment have to be taken into account, but also the cross- 
correlations between these different moments, because they are not independent of each 
other. In fact, if we consider the full correlation matrix for a combination of moments, for a 
range of sq values, it turns out to have zero eigenvalues at machine precision because of the 
strong cross-correlations. This means that standard fits, employing the full correlation 
matrix in constructing the function to be minimized, are not possible in this case. This 
remains true if we "thin out" the data, i.e., if we use fewer values of Sq on a given interval, 
by a factor two to four. This puts standard fits out of reach for simultaneous fits to 
multiple moments, forcing us to either use a different fitting strategy, or to drop such fits 
from consideration. 

Because of this problem, we will perform fits to multiple moments using a different "fit 
quality" Q^. will be a positive-definite quadratic form in the differences between data 
and theory; for a description of some possible choices, see Sec. IV Bl Any such can be 
minimized to give an estimate of the fit parameters, as long as we have a reasonable way 
to estimate the parameter errors and covariances associated with the fit. In this article, we 
will estimate the parameter error matrix by propagating errors through a linear fluctuation 
analysis starting from the full data covariance matrix; for details, see App. jX] 



We would like to thank Sven Menke for making the data files available to us. 

These numbers are at the right edges of the bins centered at s = 1.488 GeV^ and s ~ 3.120 GeV^, 
respectively. 
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B. Fitting strategies 



In this subsection, we explain three different fitting strategies we have used, in view of 
the problem of strong correlations described in Sec. IV A[ 

1. The simplest fit we can perform is a standard fit to a single moment. We will choose 
our single-moment fit to be the one with weight function Wq{x) = 1, since it does not 
suppress contributions from any part of the spectrum, and is sensitive enough to the 
DV part of our fitting function to give reasonably good fits for the DV parameters. 
It should be noted here that our main goal is to minimize the fit error on as- The 
only reason one cares about the "nuisance" parameters kv,a, Jva^ ctv,A, and f3v,A is 
that they describe part of the physics, and as such have to be taken into account in 
any fit. We emphasize again that, under the assumption that our DV ansatz ( I2.15P 
gives a good description of the DVs, there are no reasons to limit ourselves to pinched 
weights and, in fact, strong arguments not to do so. 

We have also considered fitting the spectral function directly, since it is maximally 
sensitive to DVs in the kinematically allowed region. Such a fit can be cast in terms 
of an FESR obtained by choosing w{s) = 1 and replacing p^^~^^\s) with its derivative 
on the left-hand side of Eq. f l2.3p . It turns out that it is not possible to determine 
from such a fit, basically because the spectral function is much less sensitive to than 
it is to the DV parameters. In contrast, pinched weights suppress the contribution 
from DVs more than wq does. We find, in fact, that fits with a single doubly-pinched 
weight are not stable if one tries to fit both the OPE and (strongly suppressed) DV 
parameters. In our experience, the most stable and precise results from a single- 
moment fit are obtained using Wq{x) = 1. Fits to singly-pinched weights also appear 
to work well: we have checked standard fits with weights ^2(2^) or 1 — and find 
results in excellent agreement with those reported in Tables [T] and S] below, though 
with somewhat larger errors on a^. 

2. It is of course interesting to see whether simultaneous fits to multiple moments can be 
used to reduce errors, in particular the error on ag- However, in this case, we run into 
the problem of strong correlations described above in Sec. IV A[ The simplest solution 
is to omit correlations in constructing the fit quality Q^, and choose a which is 
diagonal in the differences between data and theory. Working with a set of sq, {sq} 
in some fitting window, and letting 5I^^\so) be the error on the weighted spectral 
integral I^'^\so), obtained using the full data covariance matrix, such a fit quality has 
the form 

where we have made the dependence of the weighted theory integral on the set of fit 
parameters p explicit and the outer sum runs over the set of weights included in the 
analysis. 



■^^ In Ref. [l3| fits to the spectral function were presented, but there as was kept fixed. Note that contributions 



to pv/A from the higher D terms in the OPE are negligibly small. 
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Often, such a fit is referred to as "uncorrelated." Indeed, if would be interpreted 
as a the standard errors obtained from such a fit would miss the effect of 
correlations and be (significantly) underestimated. However, we emphasize that we 
will not compute parameter errors in this way; instead we will propagate errors using 
the linear fluctuation analysis of App. |X1 thus taking into account all correlations 
explicitly. 

While we will not report on fits to Eq. f lS.ip . but instead rely on fits described under 
items 1 and 3, we have carried out many such fits. They yield results fully consistent 
with those we do report, but always lead to larger parameter errors. 

3. A third type of fit we will consider is one in which incorporates the correlation 
sub-matrix corresponding to each individual moment employed in the fit, but not 
the cross-correlations between different moments. Full correlations are again to be 
included via the linear fluctuation analysis described in App. |X1 We choose 

QLc. = E E (4"^n^o) - 4"^(^o;p)) (c^'^O • (^"'^(^o) - 4"^(^o;p)) , (5.2) 

with the covariance matrix of the set of moments with fixed weight w and sq 
running over the chosen fit window range. The motivation for this form is that the 
cross-correlations between two moments arise mainly because the weight functions 
used in multiple-moment fits are in practice close to being linearly dependent (even 
though, as a set of polynomials, of course they are not). This dependency might be 
reinforced by the relatively large errors on the data for values of s toward m^, because 
it is primarily in this region that the weights Wq, W2 and differ from each other. 

A key observation is that it does not matter which fit quality one chooses, ^"^ as long 
as errors are propagated appropriately. Whatever the motivation for a particular choice, 
such a choice is useful if it turns out to allow a reliable fit, and to reduce errors on the 
fit parameters. We note that, of course, it is not possible to use the minimum value of 
diagonal or block diagonal of Eqs. ( 15. ip and (15. 2p obtained in such a fit in order to derive 
a confidence level; only the relative size of minimum values compared between different fits 
with the same choice of is meaningful. 

VI. FITS 

In this section, we will present the results from our fits, using the parametrization of the 
theory explained in Sec. IIVI and employing the strategies of Sec. IV B[ All fits are based 
on the original, unmodified OPAL data, including the OPAL normalization for the vr-pole 
contribution, which corresponds to a central value of 94.0 MeV for f^. 

Section IVl Al contains our "benchmark" fit, which is a standard fit of ^ex°K^o) the 
V channel. In this case, the fit quality is the standard function, which of course employs 
the full 4(7) (so) covariance matrix, generated from the covariance matrix of the original 



We find indeed that this leads to much larger errors than naive "x^" errors would suggest. 
As long as is a positive-definite quadratic form in the differences IcxHsq) — Ith\sQ',p)- 
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Smin 


dof 


xVdof 


ds 




IV 


av 


^v 


1.3 


53 


0.44 


0.320(23) 


0.026(18) 


0.42(50) 


0.54(54) 


2.85(30) 


1.4 


50 


0.35 


0.311(19) 


0.019(13) 


0.23(44) 


-0.29(64) 


3.27(33) 


1.5 


47 


0.36 


0.307(18) 


0.017(11) 


0.16(42) 


-0.52(74) 


3.38(38) 


1.6 


44 


0.38 


0.308(20) 


0.018(15) 


0.22(51) 


-0.47(82) 


3.36(41) 


1.7 


41 


0.38 


0.305(19) 


0.012(12) 


0.03(51) 


-0.61(86) 


3.41(41) 


1.3 


53 


0.45 


0.332(37) 


0.037(27) 


0.64(53) 


0.57(58) 


2.80(33) 


1.4 


50 


0.36 


0.326(27) 


0.023(16) 


0.35(48) 


-0.32(64) 


3.27(33) 


1.5 


47 


0.37 


0.322(25) 


0.020(13) 


0.25(44) 


-0.57(73) 


3.39(38) 


1.6 


44 


0.38 


0.323(27) 


0.022(20) 


0.31(57) 


-0.53(81) 


3.38(41) 


1.7 


41 


0.39 


0.320(25) 


0.014(13) 


0.08(53) 


-0.68(85) 


3.43(40) 



TABLE 1: Standard fits to Eq. / fO) with w{s) = I, V channel. FOPT results are shown 

above the double horizontal line, CIPT results below. Errors are standard errors; and /3y in 
GeV"2 

OPAL data. We will also consider a combined fit of the same moment to the V and A 
channels. 

In Sec. I VI Bl we consider simultaneous fits to I^"^^^ and again for both the 

pure V channel and combined V^A channel cases. As already discussed in Sec. IV At we 
find that standard fits are not possible. Our best results in the V channel originate from 
minimizing the fit quality (15. 2p . Our main conclusion is that, presumably because of the 
strong correlations between different moments, these type of fits do not help reduce the error 
on as significantly. They do, however, provide cross-checks, verifying that our ansatz f l2.15p 
also describes moments other than just including the non-strange component of R^, 

Rv+A;ud{so) = Rv;udiSo) + RA;udi'^o) , (6.1) 

which is proportional to I^xV+Ai^o)- They also give access to the V- and A-channel OPE 
coefficients C^y, Cq^a, Csy and Cs,a- 

In Sec. IVI CI we consider the additional errors originating from the truncation of per- 
turbation theory, and in Sec. IVI Dl we show that our fits both give a good description of 
Rv+A;ud, and satisfy, within errors, the Weinberg and DGMLY V — A sum-rule constraints. 

A. Fits with the weight wo{x) 

We begin with a standard fit to the \^-channel w{s) = wo{x) = 1 FESR. Fit results 
are presented in Table [T|, which shows all fit parameters, as well as the number of degrees 
of freedom (dof), and the per degree of freedom. Errors are standard errors; errors 
computed with Eq. flA5p are typically somewhat larger, but similar in size. 

We observe that there is excellent stability for the results with s^m = 1-4, 1.5 and 
1.6 GeV^, with the errors getting somewhat larger at 1.6 GeV^. At values of s^m > 1-7 GeV^, 
7y becomes very small and tends to go negative, which is clearly unphysical. However, within 
errors such fits are always consistent with those shown in the table. We choose the results 
obtained at Smin = 1.5 GeV^ to fix our central values, and treat the spread of the fit results 



15 



0.14 




1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

^0 (GeV2) so (GeV2) 



FIG. 2: Left panel: comparison of lex (sq) and {sq) for the Smin = 1-5 GeV^ V- channel fits of 
Table [11 Right panel: comparison of the theoretical spectral function resulting from this fit with the 
experimental results. CIPT fits are shown in red (dashed) and FOPT in blue (solid). The (much 
flatter) black curves display only the OPE parts of the FOPT (solid) and CIPT (dashed) fit results. 
The vertical dashed line indicates the location of Smin ■ 



for Smin ranging from 1.4 to 1.6 GeV^ as an error. From this simple fit, we obtain for at 
the r mass the values 

a,(m2) = 0.307 ±0.018 ±0.004 (FOPT), (6.2) 
asiml) = 0.322 ± 0.025 ± 0.004 (CIPT) , 

where the second error represents the variation with Smin discussed above. We note that 
there is good stability over the full range of Smin values covered in Table [1] Since our fitting 
function is non-linear, in general errors are expected to be asymmetric. We have therefore 
also computed asymmetric errors for all the fit parameters. We find that the error on 
is nearly symmetric, and that only errors on Ky and 7y show a significant asymmetry. For 
instance, we find, at s^m = 1-5 GeV^, that, for FOPT, 

asiml) = 0.307^S:Si« , (6.3) 
Kv = 0.017^S:|]^7 ^ 

7y = O.Wt'oil GeV-2 , 
av = -0.52 1°:^^ , 

= 3.38l°:3^° GeV-2 . 

This shows that omitting DVs from the fit, which is equivalent to setting Ky = 0, would 
lead to a poor fit. For CIPT the asymmetries show the same pattern. 

In Fig. [21 we show the quality of the match between the fitted ll^''\so) and /^^"^(sq), as 
well as of the match between the experimental spectral function and its theoretical coun- 
terpart in which parameter values are obtained from the FESR fit for Smin = 1-5 GeV^. 
The left panel shows the FESR fit itself, while the right panel shows the results of this 
comparison for the spectral function case. We emphasize that the spectral function was not 
part of the fit; the agreement between the theoretical and experimental versions of pv{s) 
is an output. Agreement with data is good in the full fit window Sq > Smin = 1-5 GeV^. 
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Us 




IV 


ay 


/3y 


Us 


1 


-0.69 


-0.67 


0.70 


-0.62 


KV 


-0.69 


1 


0.99 


-0.47 


0.43 


iv 


-0.67 


0.99 


1 


-0.48 


0.43 


ay 


0.70 


-0.47 


-0.68 


1 


-0.98 


fiv 


-0.62 


0.43 


0.43 


-0.98 


1 



TABLE 2: Parameter correlation matrix for the FOPT fit with Smin = 1-5 GeV^ shown in TableUl 



The black curves, in contrast, show the OPE parts of the theoretical curves, i.e., the curves 
obtained by removing the DV contributions from the blue and red curves. It is clear that 
DVs are needed to give a good description of the data for J^!^°) and the spectral function 
itself, in agreement with our conclusion based on Eq. (16. 3p . 

There are strong correlations between the parameters shown in Table [T] and Eq. (16.31) . 
Such strong correlations are present in all our fits, and are unavoidable, given the number of 
fit parameters. We emphasize that this cannot be resolved by simply omitting the duality- 
violating part from the theory - one cannot "improve" fits by throwing out physics that is 
known to have an impact on those fits! In Table [2] we show the full parameter correlation 
matrix corresponding to the FOPT result quoted in Eq. (16.21) . The analogous matrix for the 
CIPT result looks very similar. 

We have also considered fits like those shown in Table [1], but including the contribution 
coming from the logarithmic dependence of C4 on s. For the latter, we estimated the 



quark-condensate contribution from the Gell-Mann-Oakes-Renner relation [35|, and we took 
{O'sG^uG^'^) = 0.021 GeV^. We find that corresponding changes in the numbers in Tabled] 
are at most a tiny fraction of the fitting errors. 

We performed a similar type of fit to the combination of V and A channels. If we use 
all possible sq values, we find that the standard fit function (involving the very strongly 
correlated spectral integral covariance matrix) is very flat, admitting not just "physical" 
solutions (consistent with those in Tabled]), but also solutions that are clearly unphysical 
(with, for instance, values for as drifting down to unacceptably low values as a function of 
Smin)- This happens for both CIPT and FOPT. Parameter errors for these solutions can be 
very large, consistent with the flatness of the landscape, and there is a strong sensitivity 
of central values to initial guesses for the parameter values. 

However, if we thin out the Sq values i.e., use only every nth value of IfJ.^{so), for some 
value of n > 1, we find that fits with n = 2, 3 or 4 are much more stable than the one 
with n = 1 (no thinning) described above. In Table [3] we show our results for n = 3, which 
is the choice leading to the most stable fits.^^ The problem disappears when we use fit 
quality (15. ip . computing errors with Eq. (1A5I) . but this method leads to significantly larger 
errors. 

Choosing again the fit with s^m = 1-5 GeV^, we find for our combined V and A channel 



18 Negative values for jv.a can appear because we use the analytic form of the integral in Eq. (12.141) . Note, 
however, that all values we find in the fits are consistent with a positive value, as physically required. 
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Smin 


dof 


xVdof 






IV A 




/3v,A 


1.3 


30 


0.81 


0.326(13) 


0.0186(88) 
0.094(51) 


0.18(35) 
1.15(35) 


0.35(46) 
0.21(79) 


2.95(27) 
-3.42(45) 


1.4 


28 


0.69 


0.317(15) 


0.0140(68) 
0.085(39) 


0.01(33) 
1.06(30) 


-0.31(58) 
-0.5(1.1) 


3.29(31) 
-3.06(61) 


1.5 


26 


0.69 


0.308(16) 


0.0134(69) 
0.110(71) 


-0.01(34) 
1.15(35) 


-0.67(70) 
-1.1(1.1) 


3.46(36) 
-2.73(59) 


1.6 


24 


0.73 


0.308(17) 


0.0150(98) 
0.15(13) 


0.06(41) 
1.29(45) 


-0.64(74) 
-1.2(1.2) 


3.45(38) 
-2.67(65) 


1.7 


22 


0.68 


0.304(18) 


0.0131(99) 
1.0(2.0) 


0.00(42) 
2.14(84) 


-0.80(77) 
-2.3(2.1) 


3.51(38) 
-2.1(1.1) 


1.3 


30 


0.85 


0.346(19) 


0.026(12) 
0.072(35) 


0.37(34) 
1.01(31) 


0.40(51) 
0.10(85) 


2.89(30) 
-3.38(48) 


1.4 


28 


0.72 


0.336(21) 


0.0170(86) 
0.075(33) 


0.10(35) 
0.99(28) 


-0.36(56) 
-0.5(1.1) 


3.30(30) 
-3.04(61) 


1.5 


26 


0.71 


0.325(22) 


0.0152(80) 
0.101(66) 


0.05(35) 
1.11(35) 


-0.73(67) 
-1.2(1.0) 


3.48(35) 
-2.74(58) 


1.6 


24 


0.75 


0.324(23) 


0.017(12) 
0.14(12) 


0.11(44) 
1.25(45) 


-0.72(71) 
-1.3(1.1) 


3.47(37) 
-2.66(62) 


1.7 


22 


0.70 


0.318(23) 


0.014(11) 
1.0(1.9) 


0.03(44) 
2.11(86) 


-0.88(75) 
-2.3(2.1) 


3.54(37) 
-2.1(1.1) 



TABLE 3: Standard fits to Eq. \2. 3\) for w{s) = 1, combined V&A channels. FOPT results are 
shown above the double horizontal line, CIPT results below. The first line for each Smm gives the 
V DV parameters; the second line the A ones. Every third value of sq starting at Smin is included 
in the fits. Errors are standard errors; jv,A CLn-d fSy^A in GeV~^. 

fit the results 

as(m2) = 0.308 ±0.016 ±0.009 (FOPT), (6.4) 
as{ml) = 0.325 ±0.022 ±0.011 (CIPT) , 

with the second error representing, as above, the variation over the neighboring Smin values, 
1.4 and 1.6 GeV^. The results reported in Tables [T] and [3] are in good agreement. We note 
that the central values for jy shown in Table [3] are very small, compared with those in 
Table [H but given the errors there is no inconsistency. In Fig. [3] we show the quality of the 
fit in the panels on the left, for V (top) and A (bottom), for Smin = 1.5 GeV^. In the panels 
on the right we show again the match between the experimental spectral functions and their 
theoretical counterparts with parameter values obtained from the Smin = 1.5 GeV^ FESR 
fit. As before, black curves show only the OPE parts of the theoretical curves. Again, we see 
that the data clearly confirm the presence of DVs, which are well described by our ansatz. 
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FIG. 3: Left panels: comparison of Iex°\so) and lj:^°\sQ) for the Smin = 1-5 GeV^ combined 
A fits of Tahle\^(top: V channel, bottom: A channel). Right panels: comparison of the theoretical 
spectral function resulting from this fit with the experimental results. CIPT fits are shown in red 
(dashed) and FOPT in blue (solid). The (much flatter) black curves display only the OPE parts of 
the FOPT (solid) and CIPT (dashed) fit results. The vertical dashed lines indicate the location of 



B. Multiple-weight fits with the weights wq, W2 and 

One would naively expect that by using more moments, more information could be ex- 
tracted from the data. This would help reducing the errors reported in Tables [T] and ID 
Higher-degree weights, however, also require the introduction of additional OPE fit param- 
eters. This, in combination with the very strong correlations, may turn out to reduce the 
extra constraints placed on the parameters (ctg and the DV parameters) entering the Wq 
FESR by the additional moments. In addition, as we already pointed out in Sec. IV At it 
appears to be impossible to perform standard fits to multiple moments. 

Table m shows the results of simultaneous fits to moments with weights wq, W2, and 
for the V channel using the fit quality (15. 2p . with errors computed from Eq. flA5|) : see Fig. H] 
for a visual representation of the qualities of the resulting fits for the case Smin = 1-5 GeV^. 
Results are consistent with those presented in Sec. IVI A[ but errors are slightly larger. 

Again, the black curves in Fig. Hlshow the OPE parts of the theoretical curves, i.e., the 
results obtained by removing the DV contributions from the blue and red curve results. For 
the spectral function, 1^'^°^ and it is again clear that no good description of the data can 
be obtained without a model for DVs. This is not the case for the doubly-pinched moment 
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TABLE 4: Fits to Eq. i2.3\) with weights t/)o,2,3; V channel, using fit quality \5.S\) . FOPT results 
are shown above the double horizontal line, CIPT fits below. Errors have been computed using 
Eq. (X^) : -iy and (3v in GeV"^, Cey in GeV^ and Csy in GeV^. 



/^fs). In this case, one would expect that a reasonably good fit can be obtained without 
DVs, for values of s^m down to somewhere below ~ 2 GeV^. This is consistent with the 
results of Ref. T2| for various doubly pinched weights, and, for the doubly pinched kinematic 



weight, also with the results of Refs. js], lof; in those cases, reasonably good matches for the 
sum of the vector and axial channels were obtained without the inclusion of DVs. However, 
one would also expect that a best fit can only be obtained by shifting the OPE parameters 
relative to those reported in Table |H 

An estimate similar to Eq. fl6.2p . using the fit with Smin = 1-5 GeV^, leads to 

a,(m2) = 0.304 ±0.019 ±0.001 (FOPT), (6.5) 
asiml) = 0.322 ± 0.031 ± 0.005 (CIPT) , 

which is in excellent agreement with Eq. (16.21) . We see, however, that the simultaneous fit 
to multiple moments does not help reduce the error. It is an interesting question whether 
fit qualities other than those of Eqs. (15. ip and (15.20 exist that would lead to smaller errors. 
Fits like those reported in Table H] using fit quality (15.10 lead to consistent results, but with 
significantly larger errors than those shown in the table. 

Finally, in Table [5l we show the results of combined VSzA channel fits, similar to those 
of Table [3], but for the weights wq, W2 and w^. From this table, we obtain 

a,(m2) = 0.302 ±0.015 ±0.001 (FOPT), (6.6) 
asiml) = 0.322 ±0.024 ±0.008 (CIPT) , 

with the second error again reflecting the variation over the range Smm = 1-4 to 1.6 GeV^. 
Comparing Table with Table H] we see that adding the A channel gives little extra infor- 
mation (apart from estimates of the axial OPE and DV parameters). Adding the A channel 
increases the central values of jy, but all parameter values are consistent between these two 
tables within (sometimes substantial) errors. The fit with Smin = 1-7 GeV^ is clearly not 
meaningful (in particular for FOPT), and the errors indicate that at this value of Smin, 
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FIG. 4: V -channel fits of Table^ showing the theoretical and experimental versions of the moments 
/('^o) (top left), I^'^^) (^tgp rigf^tj andl^'^''^ (bottom left), for Smin = 1-5 GeV^ Bottom-right panel- 
comparison of the theoretical spectral function resulting from this fit with the experimental results. 
CIPT fits are shown in red (dashed) and FOPT in blue (solid). The (flatter) black curves display 
only the OPE parts of the FOPT (solid) and CIPT (dashed) fit results. The vertical dashed lines 
indicate the location of Smin ■ 



is very flat in some directions in parameter space. Indeed, restricting the value of as to 
the value obtained at Smin = 1-6 GeV^ leads to a good flt for the remaining parameters, 
consistent with the unrestricted s^m = 1-6 GeV^ flt, and with a value for almost equal 
to the value reported in the table. 

Let us next deduce some implications of the above flts for the breaking of the factorization 
hypothesis in the D = 6 condensates. To begin with, we will assume that the D = 6 conden- 
sates are dominated by their leading-order contribution. The corresponding contribution to 
the V/A correlators is given by |2|: 



Ce,v/A = -8 vr^a, ((^7^ {^^j t^d) (d^ {^^ t^u)) (6.7) 
E {{ui^t''u + d^^t''d){qrty)) + 0{al) , 

q=u,d,s 

where, in the flrst line, the upper Dirac structure 75 corresponds to the V channel and the 
lower 1 to the A channel. The two four-quark condensates can be parametrized in terms of 
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TABLE 5: Fits to Eq. ^2. 3\) with weights wq^2,3, combined V and A channels, using fit quality 115. 
FOPT results are shown above the double horizontal line, CIPT fits below. Errors have been 
computed with Eq. M5|) ." 7y,A md fiy^A in GeV~^, Cey and Cq^a in- GeV^ and Cgy and Cg,^ in 
GeV*. The first line for each Smin gives the V channel OPE and DV parameters; the second line 
the A channel ones. Every third value of sq starting at Smin is included in the fits. 

their factorization values by 

4 _ 

{{qaf.t''qj){qjYt''qi)) = -- {qiqi){qjqj)- pi , (6.8) 

4 _ 

{{qai^i5t''qj){qji''i5t''qi)) = g{qiqi){qjqj)-p5 , 

where the parameters pi and ps would be equal to one if the vacuum-saturation approxima- 
tion were exact. Further assuming that isospin breaking is small in the light u- and d-quark 
sector, that is (uu) = {dd) = (qq), Eqs. (16. 7p and 16.81 imply 

C^6yM = |vr^a.(gg)^('^;[;'^^). (6.9) 

Inverting Eq. (16. 9p . on the basis of the results of Tabled estimates of the parameters pi,5 
can be deduced. As representative examples, for the central fits with Smin = 1-5 GeV^, one 
obtains: 

pi = -1.4±3.2, p5 = 3.3±1.5 (FOPT), (6.10) 
pi = - 0.9 ± 3.1 , p5 = 3.4 ± 1.5 (CIPT) , 
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where {qq){ml) = — (272MeV)^ |33|, together with our results for as{m'^), has been 
employed. The central results of Eq. (16.101) display sizable deviations from the factorization 
values pi = p5 = 1, though, given the large uncertainties, the significance is not very high. 
The employed perturbative resummation scheme does not seem to play a big role for the 
condensate or DV parameters, suggesting that this choice is mostly compensated for by the 
differing values. 

C. Errors from truncating perturbation theory 

One of the uncertainties afflicting any determination of as is that perturbation theory 
needs to be truncated, irrespective of whether the CIPT or FOPT resummation scheme is 
used for the truncated series. As already mentioned in Sec. [ITl we use the known values of 
Cni up to n = 4, together with the estimate C51 = 283. We assign a 100% error, ±283, to 
this estimate, using this as a measure of the truncation uncertainty. 

For the fits presented in Eqs. f l6.2p . f l6.4p . fl6.5p and f l6.6p we find that this variation of C51 
leads to a shift of at most ±0.006 in Q;s(m^); for Eq. (16. 2p the shift is ±0.005. The shifts in 
other parameters are also small (well within fitting errors). 



D. Consistency with the V + A and V — A chiral sum rules 



While Figs. |2H1] show that the parameter values obtained from the fits corresponding to 
these figures give a good description of the data, with those parameter values in hand one 
may also consider other quantities. The total non-strange scaled V + A branching fraction 
Rv+A;ud {cf. Eqs. (12. 5p . (12. 7p and (16. ip ) has always played a central role in the study of 
hadronic r decays. In particular, following Refs. [sllsl, we may consider Rv+A;ud{so) for 
a hypothetical r of mass-squared = sq, as a function of Sq. We show our version of 
this quantity in Fig. El using the parameter values for the fit with Smin = 1.5 GeV^ of 
Table El the only fit reported that simultaneously yields all parameters needed to evaluate 
Rv+A;udi^o)-^^ Clearly, our result compares well with the fits shown in Fig. 10 of Ref. [6[,^° 
especially keeping in mind that in Fig. [5] we only show errors on the experimental spectral 
integrals, and not on the theory curves. Both CIPT and FOPT describe the data well down 
to So = 1.5 GeV^. We used the same values for Sew = 1.0194 and iKidp = 0.9512 as 
Ref. (sf in order to plot Rv+A;ud{so). In view of the discussion in Ref. [12] of the analysis 
of Refs. jBI, 0|, we conclude that our fits pass the test of Rv+A;udi^o) much better than the 
original analyses of Refs. js], |^ , and over a wider range of Sq than the alternate, self-consistent 
fits obtained ignoring DVs in Ref. [l2j. We emphasize, though, that it is not sufficient to 
find a satisfactory description of this quantity only - at the very least all FESRs used in the 
fits should show a similarly good match between experiment and theory as a function of Sq, 
as was shown to be the case for our fits in Sees. IVI Al and IVI Bl above. 

We may perform a similar test on the fits reported in Table HI Of course, in that case 
only l^-channel parameters are available, so one should consider the corresponding ratio 
Rv;ud as a function of sq. For the V channel, Rv;ud{^o) coincides (up to the overall factor 



We included the pion-pole contribution to the longitudinal part in Eq. p.7p in both the data points and 
the theory curves in Fig. [5] 



For an early investigation of this type, see Ref. 37| . 
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FIG. 5: Rv+A;udiso) OS a function of sq, with Smin = 1-5 GeV^ theory curves from Table\^ for 
CIPT (dashed red curve) and FOPT (solid blue curve). 
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FIG. 6: First (left, in GeV^j and second (right, in GeV^J Weinberg sum rules. Data are used for 
s < Ssw, while the DV ansatz 1^2. 15\) with values from Table\^for Smin = 1-5 GeV^ has been used 
for s > Ssw FOPT figures are shown; CIPT figures look essentially identical. 



127r'^ SEw\Vud\'^) with /^|^^)(so). The corresponding match between data and theory for Smin = 
1.5 GeV^ is shown in the left bottom panel of Fig. |H 

We can also test our results by considering how well they satisfy the classical chiral V — A 
sum rule constraints represented by the two Weinberg sum rules [l5[ and the DGMLY sum 
rule for the vr electromagnetic mass splitting [l^. These tests focus specifically on the DV 
part of our spectral functions, since the OPE for the V — A correlator (and hence the OPE 
part of the spectral ansatz) has no = contribution and D > 2 OPE contributions to 
Pv/a{s) are tiny, and can be ignored. 

Weinberg's sum rules can be written as 

f ds {pr'\s) - pT\s)) = [ds[p^i\s)-p^^{s))-2fl = Q, (6.11) 
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where we assumed, as before, that we can neglect terms of order miirij with i,j = u, d, even 
though there is a term of order rriimja'^ hnearly divergent in Sq in the second sum rule. The 
fact that this term is still very small at Sq = amounts to the observation that the chiral 
symmetry breaking terms in the second of Eq. (16. lip are not visible in the data. This means 
that in our test of the second sum rule we can assume ourselves to be effectively in the chiral 
limit, in which the divergence does not appear. In Fig. [6] we show the integrals in Eq. (16. lip 
as a function of the "switch point" Ssw below which experimental spectral data is used and 
above which the difference of the V and A DV ansdtze of Eq. (I2.15P is employed for pv-a{s). 
The DV contributions were obtained using the DV parameter values of Table |5l If the DV 
ansatz is, as we have assumed, reliable in the window of sq employed in the FESR fits which 
produce these DV parameter values, the V — A sum rules should be satisfied for all values 
of Ssw lying in this sq fit window. We see that this condition is well satisfied for both of the 
Weinberg sum rules. The errors shown are those from the experimental {s < Sgw) pai't of 
the integral on the left-hand side of the sum rules only. 

The first Weinberg sum rule is, of course, closely related to the difference of the V and A 
Wo{x) = 1 FESRs. Therefore, the very good quality of the matches between the experimental 
spectral integrals and our theoretical representations thereof precludes the first Weinberg 
sum rule being badly broken by our fits. The second Weinberg sum rule (as well as the 
DGMLY sum rule discussed below) can be viewed as a prediction, because the fits of Table [5] 
do not involve any weight with a term linear in s. We also note that if indeed we take the 
second sum rule in the chiral limit, we should omit the term — 2m^/^. The value of this 
term is —0.00034 GeV^, and it can thus indeed safely be dropped from the sum rule — the 
difference would not be visible in the figure. 

Finally, we consider the DGMLY sum rule for the vr electromagnetic mass difference. To 
leading order in the chiral expansion, one has that 

l^dss\og{s/p'){p^^\s)-p^}\s)) = -^{ml,-mlo) , (6.12) 

where a is the fine-structure constant, and /o the vr decay constant in the two- flavor chiral 
limit. On the right-hand side we take as input the values /o = 87.0 ± 0.6 MeV 38|] and 
1711+ — rriio = 0.00126 ± 0.00008 GeV^ |l9|.^^ Because of the second Weinberg sum rule. 



the left-hand side of Eq. (I6.12p is, in fact, independent of the scale p. In Fig. [7] we show 
the left-hand and the right-hand sides of Eq. (16.120 . as a function of Sgw The left-hand 
side is represented by the data points, while the gray band represents the right-hand side, 
including the error resulting from the uncertainties in /o and the pion electromagnetic mass 
difference. As for the Weinberg sum rules, the integral on the left-hand side is computed 
using experimental data for s < Sgw, and the DV ansatz, with DV parameters from the 
Smin = 1-5 GeV^ entries of Table |5l for s > Ssw The errors shown again come from the 
experimental data part of the integral only. We chose p^ = 2.5 GeV^, which makes the 
errors in the figure relatively small. For lower values of p"^ the errors are larger for larger 
Ssw, but the results are consistent with Eq. (I6.12p being satisfied for all p'^ between 1.5 GeV^ 
and ml. 



Our uncertainty on the electromagnetic contribution to the difference to^+ — m^o comes from an estimate 



of the contribution of m„ — md to the pion mass difference 39 1 
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FIG. 7: The DGMLY vr electromagnetic mass difference sum rule: left-hand side (data points) 
and right-hand side (gray band) of Eq. h6.1S\) . as a function of Ssw, in GeV^. Data are used for 
s < Ssw, the DV ansatz ^2.15\) . with fit parameter values from the Smin = 1-5 GeV^ entries of 
Tahle\^ for s > Sgyj. The FOPT figure is shown; the CIPT figure is essentially identical. The gray 
band represents an estimate of the error on the right-hand side of Eq. h6.10) . 



One might consider using these sum rules as a further constraint on the DV parameters. 
Of course, this can only be done for combined VhA fits, which requires us to assume that 
the larger of the two Smin values for the V and A channels still lies sufficiently far below 
to make such a combined fit, using our DV ansatz in both channels, reliable. Since in this 
work our primary results are obtained from purely \^-channel fits, which do not require this 
assumption, we postpone this possibility to a future investigation. 



VII. SUMMARY OF RESULTS 

Our most stable results come from fits to only the V channel. Moreover, since we need 
to include a model for DVs, using only the V channel avoids the additional assumption that 
our ansate already adequately describes the A-channel data in the interval between 1.5 GeV^ 
and m^. We emphasize, however, that all the results from our combined V and A channel 
fits are consistent with those from our V^-channel analyses. In addition, they pass several 
further V + A and V — A channel tests, as demonstrated in Sec. IVIDI 

In view of these observations, we will choose the fit to in the vector channel as our 
central result. Adding the errors from the fit, the variation of s^m and the variation of C51 



See, for instance, Ref. [36|, where, however, the difference V — A was modeled with an expression of the 
form (|2.15p . Given the results of Tables [3] and [5l such a parametrization does not seem to be favored by 
the data. In the present work we avoid any additional assumptions about the relations between the DV 
parameters in the V and A channels by introducing separate DV parameters for each channel. 
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in quadrature, we find 



asiml) = 0.307 ±0.019 {MS, nj = 3, FOPT) , (7.1) 
asiml) = 0.322 ± 0.026 (MS, nj = 3, CIPT) . 

We do not average over CIPT and FOPT results, because we believe that it is useful to see 
the difference between values for obtained with the two resummation schemes after all 
non-perturbative effects have been consistently taken into account. Running these values 
up to the Z mass Mz yields (iot^^ 

a,(M|) = 0.1169 ±0.0025 (MS, nj = 5, FOPT) , (7.2) 
a^(M|) = 0.1187 ±0.0032 (MS, n/ = 5, CIPT) , 

where we symmetrized the slightly asymmetric errors one obtains after running up to the Z 
mass. 

These values can be compared to those obtained by OPAL from the same data [6]. The 
OPAL values at the r mass are 

asiml) = 0.324 ±0.014 (MS, nj = 3, FOPT, OPAL) , (7.3) 
asiml) = 0.348 ± 0.021 iWS, nj = 3, CIPT, OPAL) , 

where we added the experimental and theoretical errors quoted by OPAL in quadrature. We 
observe the shift to lower central values, with somewhat larger errors that follow from using 
our new framework for analyzing the data. We also note that our CIPT and FOPT values 
are somewhat closer, and that, because of the larger errors, the difference between our two 
values for as is less significant. 



VIII. CONCLUSION 

In this article, we provided a new framework for the extraction of and other OPE 
parameters from hadronic r decays. This new framework combines two elements that have 
not been taken into account in the "traditional" analysis of hadronic r decays. One is a 
consistent treatment of the OPE, as discussed in Ref. [l2|], and the other is a detailed quan- 
titative estimate of the effect of violations of quark-hadron duality, using a parametrization 



proposed in Refs. |27|, |29| . As explained in detail in Sees. [TTTl and HVl these two elements are 
intricately intertwined. 

Our new framework comes with the price of introducing four new fit parameters for each 
channel, the vector and axial-vector DV parameters. Nevertheless, we demonstrated that 
our method is feasible by presenting a rather complete analysis based on the OPAL data 
for the V and A spectral functions j6|. With the larger number of parameters, and the 
corresponding need to vary sq away from m^, it should come as no surprise that our errors 
are typically larger than those found in earlier analyses, which simply did not take DVs into 
account quantitatively. We emphasize that this means that the systematic errors of those 
earlier analyses were understimated - we believe significantly so in some cases. 



23 



The specifics of tlie evolution to the Z mass are as discussed in Ref. 



27 



Our analysis leads us to a new estimate for both the central value of and the error, 
which, in our opinion, should be interpreted as superseding previous estimates in the liter- 
ature; for our result, based on OPAL data, see Sec. IVIII We found our most reliable fits 
to be those of the V channel, although fits including also the A channel lead to results 
consistent with our most precise V^-channel fits. As our primary concern in this article 
is with previously underestimated non-perturbative effects, we presented results for both 
contour-improved as well as fixed-order perturbation theory. Our analysis was based on the 
original OPAL data, unmodified for subsequent changes in the various exclusive branching 
fractions. This choice was made in order to facilitate interpretation of the differences in 
our results from those obtained by OPAL. With this choice, these differences are solely the 
result of differences in the analysis method. We plan to present an analysis of OPAL data 
with updated normalizations in the near future. 

The accuracy of the results presented here depends in part on our ability to correctly 
model the physics present in duality violations. Conservatively, our results can be seen as 
providing a lower bound on the error introduced by ignoring duality violations. However, 
the stability of a^, a purely perturbative parameter, across the range of moments analyzed 
here provides strong support for the validity of our ansatz. We therefore surmise that not 
only is the ansatz able to accurately describe the data, but also that it provides a reasonable 
quantitative description of the physics of duality violations in the light-quark V and A 
channels. 

In cases with multiple weights, standard fifs were not possible, and we performed 
alternate fits, propagating errors as described in App. |X1 We point out, however, that our 
final result is based on a standard fif fo the moment All other fits yield results 

completely consistent with this result, including the error on The error on obtained 
from our 

fit to /(fo) is very close to that obtained with the method of App. |X1 Because 
our parameter covariance matrix obtained with that method scales linearly with the data 
covariance matrix, the error on ag will be reduced once the improved spectral function data 
expected from BaBar and/or Belle becomes available. 

We observe that a difference remains between the central values for obtained using 
CIPT and FOPT, though this is less significant than the difference found previously by 
OPAL, cf. Eq. (17.31) . In this context, we note that, as mentioned already in Sec. IIV Al a 
term linear in s in any of the weights employed in Eq. fl2.3p picks out the D = 4 term in the 
OPE, which parametrizes the leading renormalon ambiguity in the perturbative expansion. 
This, then, raises the question whether differences between the behavior of CIPT and FOPT 
fits, including those associated with any dependence on the choice of weight, might be used 
to constrain renormalon pole models that have been used previously to investigate the 
resummation of theperturbative series for the various sum rules employed in the study of 
hadronic r decays [sl, We plan to pursue such an investigation in a future work. 
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Appendix A: Error propagation 

Consider a fit quality 

Q'^[di-U{p)]Cojj[dj-tM , 



(Al) 



in which di are the binned data, ti{p) is a function that describes this data set for a set of 
parameters p, and Co is a positive-definite, symmetric, but otherwise arbitrary matrix. In 
this Appendix, we use the summation convention for repeated indices. The parameters p 
are determined by finding the global minimum of Q^, which satisfies 



^--2^C-.K.-t,(p-)]^0 



(A2) 



dpa dpa 

Varying, in this equation, the data by an amount 5di, and the parameters by 5pa leads to 



dpadp/s 



dpa 



0,ij 



Sdj - 



. 



(A3) 



If the fit is good, so that the deviations dj — tj{p) are small, we may ignore the term with 
the second derivative, leading to 



dpp 



(A4) 



or, for the covariance matrix (SpaSp^), 



in which 



and 



dpa 



'''' dp, ' 



Cke = {5dk5(k) 



(A5) 



(A6) 



(A7) 



is the data covariance matrix. This provides us with an estimate for the full correlation 
matrix for the parameter set p. We note that, if Co,ij is chosen to be equal to the data 
covariance matrix Cjj, this expression simplifies to 



{Spjpi3) = A J . 



(A8) 
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This is equal to the usual error matrix estimate, given by the inverse of one-half times 
the second derivative of at its minimum, if, again, the fit is good enough to ignore terms 
proportional to di — Uip). 
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